Exactly solvable pairing Hamiltonian for heavy nuclei 
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We present a new exactly solvable Hamiltonian with a separable pairing interaction and non- 
degenerate single-particle energies. It is derived from the hyperbolic family of Richardson- Gaudin 
models and possesses two free parameters, one related to an interaction cutoff and the other to the 
pairing strength. These two parameters can be adjusted to give an excellent reproduction of Gogny 
self-consistent mean-field calculations in the canonical basis. 
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PACS numbers: 02.30.Ik, 21.60.Fw, 21.60.Jz 

Pairing is one of the most important ingredients of 
the effective nuclear interaction in atomic nuclei as it 
was early recognized by Bohr, Mottelson, and Pines Q 
in an attempt to explain the large gaps observed in 
even-even nuclei. They suggested that the recently pro- 
posed Bardeen-Cooper-Schriefer (BCS) [H theory of su- 
perconductivity could be a useful tool in nuclear struc- 
ture although care should be taken with the violation 
of particle number in finite nuclei. Since then, BCS or 
the more general Hartree-Fock-Bogoliubov (HFB) the- 
ory combined with effective or phenomenological nuclear 
forces has been the standard tool to describe the low en- 
ergy properties of heavy nuclei. Improvements over BCS 
or HFB came through the restoration of broken symme- 
tries, specially particle number projection which is still 
a problem not satisfactory solved with density depen- 
dent forces 0]. From a different perspective, Richardson 
found an exact solution of the constant pairing problem 
with non-degenerate single particle energies as early as in 
1963 Though highly schematic, the constant pairing 
force has been used for decades in nuclear structure with 
several approximations (BCS, RPA, PBCS, etc.) but 
scarcely resorting to the exact solution. Almost forgot- 
ten, the exact Richardson solution was recovered within 
the framework of ultrasmall superconducting grains Q, 
in which not only number projection but also pairing fluc- 
tuations were essential to describe the disappearance of 
superconductivity as a function of the grain size. 

By combining the Richardson exact solution with the 
integrable model proposed by Gaudin Q for quantum 
spin systems, it was possible to derive three families of 
integrable models called Richardson-Gaudin (RG) mod- 
els [7[. The rational family, extensively used since then, 
contains the Richardson model as a particular exactly 
solvable Hamiltonian as well as many other exactly solv- 
able Hamiltonians of relevance in quantum optics, cold 
atom physics, quantum dots, etc. Q. However, the other 
families did not find a physical realization up to very 
recently when it was shown that the hyperbolic family 
could model a o-wave pairing Hamiltonian in a 2 dimen- 
sional lattice [of , such that it was possible to study with 



the exact solution an exotic phase diagram having a non- 
trivial top olog ical phase and a third order quantum phase 
transition [10[ . In this letter we will show that the hyper- 
bolic family give rise to a separable pairing Hamiltonian 
with 2 free parameters that can be adjusted to reproduce 
the properties of heavy nuclei as described by a Gogny 
HFB treatment. 

Let us start our derivation with the integrals of motion 
of the hyperbolic RG model Q, which can be written in 
a compact form [llj as 
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where Sf, Sf, are the three generators of the SU{2)i 
algebra of copy i with spin representation Si such that 
{Sf) = Si{si -f 1). We assume L S'J7(2)-algebra copies, 
i = 1, . . . , L. The L operators Ri contain L free param- 
eters rji plus the strength of the quadratic term 7. The 
integrals of motion ([1]) commute among themselves and 
with the z component of the total spin S*^ = J2i=i ■ 
Therefore, they have a common basis of eigenstates which 
arc parametrized by the ansatz 
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where \v) is the vacuum of the lowering operators 
S~\v) = and the Ep (/3 = 1, • • • , M) arc the pair en- 
ergies or pairons which are determined by the condition 
that the ansatz ([2|) must satisfy the eigenvalue equations 
Ri |*A/) = r, \^m) for every i. 

In the pair representation of the SU{2) algebra, the 
generators are expressed in terms of fermion creation and 
annihilation operators = c\cX = (S*")^, = (c|ci -I- 

c|c-— l)/2. Each SU{2) copy is associated with a single 
particle level i, with i the time reversed partner, and M 
is the number of active pairs. The vacuum \v) is defined 
by a set of seniorities, \v) = \vi^V2, . ■ ■ ,vi) ^ where the 



2 



seniority Vi = 0,1 is the number of unpaired particles in 
level i, which determines the spin associated to the level 
as Si = (1 — ^'i)/2. The blocking effect of the unpaired 
particles reduces the number of active levels to Lc = 

L-T,^^^■ 

Although any function of the integrals of motion 
generates an exactly solvable Hamiltonian, we will re- 
strict ourselves in this presentation to the simple lin- 
ear combination H = ^'^iJliRi- Defining A = 
(1 + 27(1 — M) + jLc)^^, and after some algebraic ma- 
nipulations the Hamiltonian reduces to 
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where G = 2A7 is a free parameter. 

This Hamiltonian, expressed in a 2 dimensional mo- 
mentum space basis gave rise to the celebrated px + ipy 
model of p-wave pairing [l^l- However, if we interpret 
the parameters rji as single particle energies correspond- 
ing to a nuclear mean-field potential, the pairing interac- 
tion has the unphysical behavior of increasing in strength 
with energy. In order to reverse this unwanted effect we 
define rji = 2{ei — a), where the free parameter a plays 
the role of an energy cutoff and Si is the single particle 
energy of the mean- field level i. Making use of the pair 
representation of the S'C/(2), the exactly solvable pairing 
Hamiltonian ([3]) takes the form 



H = Y.'^{cU + ^,o-^| (4) 

i 

-2G^ ^(a - Ej) (a - £i')c\cicj-Ci' , 
a' 

with eigenvectors given by ([2]) and eigenvalues 

E = 2aM + Y£,v,+YEf}. (5) 

Here the pairons correspond to a solution of the set 
of non-linear Richardson equations 



E 



r]i - Ep 



E 



E. 



■13' 



E, 



13 



Q_ 



(6) 



where Q = — ^ + M — 1. Each particular solution of 
Eq. ^ defines a unique eigenstate 

In order to get an insight into the solutions of ([5]) we 
show in figure [1] the ground-state pairon dependence on 
the pairing strength G for a schematic system of M = 10 
pairs moving in a set of L = 24 equally spaced single- 
particle levels (ci = i) and a cutoff a = 24. For G — > 
the pairons are all real and stay close to a set of M pa- 
rameters r]i (the M lowest rj' s for the G.S. configuration) 
in order to cancel the divergence in the r.h.s. of ([5]). As G 
increases the pairons move down in energy till they reach 
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FIG. 1: Real and imaginary parts of the ground-state pairons 
as a function of pairing strength, for a set of 24 equally spaced 
single-particle levels (e^ = i), a cutoff a = 24 and M = 10 
pairs. The inset shows the pairon distribution in the complex 
plane for two different pairing strengths. 



a critical value of G « 0.012 for which the two pairons 
closest to the Fermi level collapse to ~ —30. Immedi- 
ately after they acquire an imaginary part and expand 
in the complex plane as a complex-conjugate pair. The 
same phenomenon happens to the other pairons as G is 
further increased forming an arc in the complex plane as 
can be seen in the inset of Fig. [TJ Even though the be- 
havior of the pairons resembles that of the rational model 
, there are qualitative differences associated to the non- 
constant form of the pairing interaction that will turn out 
to be essential for the description of heavy nuclei. 

In what follows we will derive the two free parame- 
ters G and a of the integrable Hamiltonian (j4]) by fit- 
ting its BCS wavefuntion to a Gogny HFB calculation 
in the canonical basis. The HFB calculations with the 
Gogny force have been carried out with the standard DIS 
and the canonical basis obtained 
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parametrization 
by diagonalizing the Hartree-Fock (HE) field. The pairing 
tensor is not exactly diagonal, but we have checked that 
the off diagonal contributions are much smaller than the 
diagonal ones. In this approximation HFB in the canon- 
ical basis is equivalent to BCS. 

Due to the separable character of the integrable pair- 
ing interaction the state dependent gaps and the pairing 
tensor in the BCS approximation are 



Aj = 2GVq; - £t \/a - ei'Ui'Vi> = A^a^^, (7) 
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FIG. 2: State dependent gaps Ai, and pairing tensor UiVi 
for protons in ^-^^U and ^**Sm. Open circles are Gogny HFB 
calculations in the canonical basis while the continuous lines 
are the BCS results of the integrable Hamiltonian. 
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Note that the gaps A, and the pairing tensor mvi de- 
pend on a single gap parameter A and have a square 
root dependence on the single particle energy. Hence, 
the model has a highly restricted form for both magni- 
tudes that we will test against the Gogny gaps Ap = 
Si' ^ii,17i'''J'?'"ff pairing tensor ufvf, where Vf^jj^, 
are the matrix elements of the Gogny force in the canon- 
ical basis and {u'^v'^) is the HFB eigenvector. We take 
the single particle energies of the integrable Hamilto- 
nian from the HF energies of the Gogny HFB calculations 
and we set up an energy cutoff of 30 MeV on top of the 
Fermi energy. Occupation probabilities above this cutoff 
are lower than 10~^ and oscillate randomly. In order to 
fit the two parameters of the model a and G and to fulfill 
the BCS equations for the chemical potential ^ and the 
gap A, we solve the following three coupled equations for 
the chemical potential ^, the gap A and the parameter 
a: 
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TABLE I: Parameter values and correlation energies for pro- 
tons in ^'"^^Sm and ^^^U 
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/x), and the quasipar- 
Eq. © is the BCS 



number equation that fixes the chemical potential [i. 
Eq. ([T0| is a fitting of the Gogny pairing tensor wf 
with respect to the gap parameters A, i.e. we minimize 

'Y^i—ip-n (^?^P ~ with respect to A. Here we se- 

lect n levels above and below the Fermi energy in order 
to enhance the quality of the fit for the most correlated 
levels. We typically choose n ^ 10. Finally, Eq. (fTTj) 
fixes the interaction cutoff a by minimizing the differ- 
ences ^ - (A^ — A^Ja — Ei) between the state depen- 
dent Gogny gaps A^ and A^, with respect to a. Once 
/i, a and A are fixed, the pairing strength is determined 
fromEq. 
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As a first step in ascertain the quality of the hyperbolic 
Hamiltonian (jlj to reproduce the superfluid features of 
heavy nuclei, we show, in Fig. 2, the state dependent 
gaps Ai and the pairing tensor mVi for protons corre- 
sponding to two heavy nuclei, ^^^Sm and ^"^^U. Following 
the fitting procedure we consider all levels below SOMeV 
above the Fermi energy and solve selfconsistently equa- 
tions (|WTT|) for the chemical potential /x, the gap param- 
eter A and the interaction cutoff a. Fig. 2 shows a re- 
markable agreement between the Gogny force and the hy- 
perbolic Hamiltonian for the pairing tensor. The Gogny 
state dependent gaps exhibit large fluctuations due to the 
details of the two-body Gogny force. However, the gen- 
eral trend of the gaps is very well described by the square 
root y/ {a — Ei) of the hyperbolic model. Although ^asy 
has 50% more proton pairs than ^^^Sm the quality of 
the mapping is excellent for both nuclei. It is interesting 
to note that the rational model, leading to the constant 
pairing exactly solvable Richardson Hamiltonian, has a 
constant gap (a horizontal line) failing completely to de- 
scribe the Gogny gaps. Table I shows the number of pairs 
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FIG. 3: Pair energies (grey circles) of the exact ground state 
solution for protons in ^-^^U and ^^^Sm. The horizontal seg- 
ments in the real axis represent the parameters rji = 2{ei — a). 



M , the number of active levels L within the energy cut- 
off, the pairing strength G, the gap parameter A and the 
correlations energies for both nuclei. 

Once we have set up the procedure to define the pa- 
rameters of the hyperbolic Hamiltonian in the BCS ap- 
proximation, we are ready to explore the exact solution. 
For a general pairing Hamiltonian the dimension of the 
Hilbcrt space is given by the Binomial B{L,M). Using 
the information of Table I, the dimensions are 1.98 x 
i54gj^ g^^^ 4 83 X 1038 for 238u^ ^ell beyond 

the limits of a large scale diagonalization. However, the 
integrability of the hyperbolic Hamiltonian allows us to 
obtain the exact solution by solving the set of M non- 
linear coupled Richardson equations (|6]). The exact cor- 
relation energy shown in Table I are in both nuclei consid- 
erable greater than the mean-field results, reflecting the 
importance of beyond mean-field quantum correlations 
and number fluctuations. The exact ground state wave- 
function is completely determined by the position of the 
M pairons in the complex plane. Fig. 3 shows the exact 
ground states for both nuclei. Considering the structure 
of the pair wavefunctions ([2]) we may argue that ^ssy 
has 4 correlated Cooper pairs, while ^^^Sm has only 2. 
Further analysis of the Cooper pair wavefunction from 
the exact solutions as was carried out in for cold 
atoms and in [l^ for nuclei within the rational model is 
straightforward but beyond the scope of this letter. 

In summary, we have presented a new exactly solvable 
Hamiltonian with separable pairing interaction and non- 



degenerate single particle energies (|4]), which arises as a 
particular linear combination of the hyperbolic integrals 
of motion ([1]) . The separable form of the pairing matrix 
elements could be derived from a novel Thomas-Fermi 
approximation for a contact interaction in a square well 
potential [l5| . We have shown that the separable Hamil- 
tonian ([4]) with 2 free parameters is able to reproduce 
qualitatively the general trend of the state dependent 
gaps as described by the Gogny force in the canonical 
basis. At the same time, it reproduces accurately the 
HFB wavefunction represented by the pairing tensor. As 
such, our exactly solvable Hamiltonian is an excellent 
benchmark for testing approximations beyond HFB in 
realistic situations for even and odd nuclei. Moreover, a 
self-consistent HF plus exact pairing approach could be 
set up along the lines of Ref. [l^ for well bound nu- 
clei. The inclusion of exact T=l proton- neutron pairing 
within this self-consistent approach is also possible [17 [ . 
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